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I. INTRODUCTION 



That the Schrodinger equation can be linked to some 
underlying stochastic process is well known since long- 
time. This idea has received along the years a number 
of different formulations: from the Feynman path inte- 
gral [H , through the Bohm-Vigier model , to the Nel- 
son stochastic mechanics |J. In all these models the 
underlying stochastic process powering the random fluc- 
tuations is a Gaussian Brownian motion, and the focus of 
the interest is the (non relativistic) Schrodinger equation 
of quantum mechanics. This particular choice is under- 
standable because on the one hand the Gaussian Brow- 
nian motion is the the most natural and widely explored 
example of Markov process available, and on the other 
hand its connection with the Schrodinger equation has 
always lent the hope of a finer understanding of quan- 
tum mysteries. 

In the framework of stochastic mechanics, however, 
this standpoint can be considerably broadened since in 
fact this theory is a model for systems more general than 
quantum mechanics: a general dynamical theory of Brow- 
nian motion that can be applied to several physical prob- 
lems (ESQ- On the other hand in recent years we 
have witnessed a considerable growth of interest in non 
Gaussian stochasticprocesses, and in particular in the 
Levy processes [1, i! 0, Ei, El- This is a field that 
was initially explored in the 30's and 40's of last cen- 



* Electronic address: lcufaro@ba.infn.Ttl 
t Electronic address: |pusterla@pd.infn.it| 



tury [HQ in, but that achieved a full blossoming of 
research only in the last decades, also as a consequence of 
the tumultuous development of computing facilities. This 
interest is witnessed by the large field of the possible ap- 
plications of these more general processes from statistical 
mechanics [7] to mathematical finance [l^, [H, ■ In the 
physical field, however, the research scope is presently 
rather confined to a particular kind of Levy processes: 
the stable processes and the corresponding fractional cal- 
culus [3, 1 while in the financial domain a vastly more 
general type of processes is at present in use. For in- 
stance the possibility of widening the perspective of the 
Schrod ing er-Brownian pair has been considered also re- 
cently [20|, but the Schrodinger equation has been gen- 
eralized only to a fractional Schrodinger equation. At 
our knowledge, instead, the association of the more gen- 
eral Levy infinitely divisible processes to the Schrodinger 
equation has already been recognized as an important 
tool only once, in a paper [2l[ principally concerned with 
the analysis of the relativistic Schrodinger equation in the 
framework of the so called Schrodinger problem solved 
by means of Bernstein processes . Here, on the other 
hand, we intend to develop the Nelson stochastic me- 
chanics as a dynamical theory of the infinitely divisible 
processes, and to widen the horizon of its applications 
even to cases different from the quantum systems Q . 

The appeal of the stable distributions is justified by the 
properties of scaling and self-similarity displayed by the 
corresponding processes, but it must also be remarked 
that these distributions show a few features that partly 
impair their usefulness as empirical models. First of all 
the non gaussian stable laws always have infinite vari- 
ance. This makes them rather suspect as a realistic tool 
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and prompts the introduction of truncated stable distri- 
butions which, however, are no longer stable. Then the 
range of the x decay rates of the probability density func- 
tions can not exceed x~^, and this too introduces a par- 
ticular rigidity in these models. On the other hand the 
more general Levy processes are generated by infinitely 
divisible laws and do not necessarily show these disturb- 
ing features, but they can be more difficult to analyze 
and to simulate [H, [H, [2^ . Beside the fact that they 
do not have natural scaling properties, the probability 
density function of their increments could be explicitly 
known only at one time instant. In fact, while their 
time evolution can always be explicitly given in terms 
of characteristic functions, their marginal densities may 
not be calculable. This is a feature, however, that they 
share with most stable processes, since the probability 
density functions of the non gaussian stable laws are ex- 
plicitly known only in precious few cases. On the other 
hand some new applications in the physical domain for 
Levy, infinitely divisible but not stable processes begin 
to emerge: in particular the statistical characteristics of 
some recent model of the collective motion in the charged 
particle accelerator beams seem to point exactly in the 
direction of some kind of Student infinitely divisible pro- 
cess [i,[2i]. 

This paper is devoted to a discussion of a generalization 
of the Schrodinger equation which takes into account the 
entire family of the Levy processes: we will propose an 
equation where the infinitesimal generator of the Brown- 
ian semigroup (the Laplacian) is substituted by the more 
general generator of a Levy semigroup. As it happens this 
will be a pseudodifferential operator (as, in particular, in 
the fractional case) , and the Levy-Khintchin formula will 
give us the opportunity to write it down in the form of 
an explicit integro-differential operator by putting in ev- 
idence its continuous (Gaussian) and its jumping (non 
Gaussian) parts. The advantages of this formulation are 
many: first of all the widening of the increment laws 
from the stable to the infinitely divisible case will offer 
the possibility of having realistic, finite variances. More- 
over both the possible presence of a Gaussian component 
in the Levy-Khintchin formula, and the wide spectrum of 
decay velocities of the increment probability densities will 
afford the possibility of having models with differences 
from the usual Brownian (and usual quantum mechani- 
cal, Schrodinger) case as small as we want. In this sense 
we could speak of small corrections to the quantum me- 
chanical, Schrodinger equation. Last but not least, there 
are examples of non stable Levy processes which are con- 
nected to a particular form of the quantum, relativistic 
Schrodinger equation: an important link that was miss- 
ing in the original Nelson model. It seems in fact that - 
as was already pointed out a few years ago [U, - we 
can only recover some kind of relativistic quantum me- 
chanics if we widen the field of the underlying stochastic 
processes at least to that of the selfdeconiposable Levy 
processes. To avoid formal complications we will confine 
our discussion to the case of processes in just one spatial 



dimension: generalizations will be straightforward. 



II. A HEURISTIC DISCUSSION 

Let us start from the non relativistic, free Schrodinger 
equation associated to its propagator or Green function 
G{x,t\y,s) (see for example |1|]) 



ihdtip{x,t) = -—dlip{x,t) 
2m 



(1) 



G{x,t\y,s) 
tp{x,t) 



1 



e 2i(t-»)B/™ (2) 



y/2Tri{t - s)h/ m 

-|-oo 

G{x,t\y,s)tp{y,s)dy (3) 



and compare it with the Fokker-Planck equation of a 
Wiener process (Brownian motion) with diffusion coef- 
ficient D, pdf (probability density function) q(x, t) and 
transition pdf p{x,t\y,s) (see for example j28| ) 



dtqix,t) = Ddlq{x,t) 
p[x,t\y,s) 



. e 4(t-s)D 



(4) 
(5) 



q{x,t) 



y/4:TT{t ~ s)D 

+ 00 

p{x,t\y,s)q{y,s)dy (6) 



It is apparent that there is a simple, formal procedure 
transforming the two structures one into the other: 
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It is well known that this is just the result of an ana- 
lytic continuation in the complex plane. There are of 
course important differences between G and p. For ex- 
ample while p and q are well behaved prf/'s, G is not 
a wave function, as can be seen also from a simple di- 
mensional argument. This simple symmetry can then be 
deceptive, but a better understanding of its true mean- 
ing can be achieved either by means of the Feynman path 
integration with a free Lagrangian of the usual quadratic 
form, or through the Madelung decomposition |29|] of I 
and its subsequent stochastic mechanical model 0, 
Our aim here is to generalize to distributions other than 
Gaussian this simple shortcut from Wiener process to 
Schrodinger equation, and to analyze its most immedi- 
ate consequences. We postpone to a subsequent paper 
a more detailed analysis in the framework of stochastic 
mechanics. 

Let us see first of all what kind of role the gaussian dis- 
tribution plays in our Wiener-Schrodinger scheme. The 
pdf and the chf (characteristic function) of a Gaussian 
law AA(0, a2) 



q{x) 



V27ra2 
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satisiy the relations 



ip{u) 



+00 



q{x) e™^ dx , 



(7) 



1 /■+°° 



du . 



Then formally the propagator ^ and the transition 
pdf ^ respectively have as c/i/'s 

where now ip{u) — e^-^"^"" = g-T?i«2/2m ^j^^ 
a Gaussian law J\f{Q,2DT), and r is a time constant 
introduced in order to have dimensionless exponents. 
From ([7]) we then have 



J-00 
-1 /•+00 

t\y, S)^— Mu)^'''^'^ ^-^ui.-y) 

In both cases the starting point is the same chf ip of the 
centered, normal law Af{0 , 2Dt). Then we consider the 
chf [(/3(u)]^*~^-'/'^ of the t~ s stationary increments of the 
Wiener process, we pass to the imaginary time variables 
(t it), and finally we get from ([8]) the Schrodinger 
propagator ([2]). The time scale r incorporated in the ini- 
tial normal law disappears in the subsequent steps when 
we generate the chf of the increments. This feature is 
common to all the stable laws and is the embodiment of 
the stable processes selfsimilarity. 

The equation ([1]) and ([4]) can now be easily deduced re- 
spectively from ^ and (O. For instance from ^ and ([5]) 
we have 



—00 27r J— 00 



and then we can write 
idtip{x, t) = 

-00 27r J— 00 



D 



+00 „^ .+00 



dy ■ 



2ti 

-Ddl r^dy^i^ /^"e-^^"'(*-^)e-»(--^)du: 
J —00 ^'^ J —00 



which - but for a factor h - is the free, non relativistic 
Schrodinger equation ([1]). 

We are interested now in reproducing these well known 
steps starting with the chf of a non Gaussian law. Take 
now an infinitely divisible - in general non Gaussian ~ law 
with chf ip(u), and let ri{u) = ln(^(u) be its Ich (logarith- 
mic characteristic). In the following we will restrict us to 



centered laws, and we will justify this choice in the subse- 
quent sections. Infinite divisibility essentially is the prop- 
erty of a chf ip which guarantees that also t^*/'^ is a legit- 
imate chf for every real t. About the infinitely divisible 
laws and their intimate relation with the Le vy p rocesses 
there is a vast literature (see for example [1, Il5| |. and for 
a short introduction [lOl). The law of the increment of 
the corresponding Levy process then is [(^(u)]*'' ^■'^^ and 
its transition pdf is ^ with our ~ possibly non Gaussian 
- infinitely divisible chf. Then, following the procedure 
previously outlined for the Wiener-Schrodinger equation, 
the wave function propagator is ^ with our new if, and 
hence from ^ the time evolution is ruled by 



■il'{x,t) = 



-00 



27r 



+ 00 



The differential equation can then be deduced as in the 
Gaussian case and is 



1 



idtij{x, t) = ri{dx)ijj{x, t) 



(10) 



+ 00 



dy ■ 



2tt 



where now ln[93(9j,)] — ri{dx) is a pseudodifjerential op- 
erator with symbol ri{u) that is defined through the use 
of the Fourier transforms [TO, [O, Hi], HI] . A pseudodif- 
ferential operator L on a suitable set of functions h(x), 
is associated to a function i{u) called the symbol of L, 
and operates in the following way: if, coherently with 
the definition of the chf of a law (see [l^l), the Fourier 
transform of a function h is defined 



h{u) 



+ 00 



(11) 



then 



(Lh){x) = — r t{u)h{u)e-''^''du. (12) 
When the symbol is a polynomial of degree n 

n 

£{u) = ^ ak{iu)'^ 

k=l 

then L is a simple differential operator of order n 

n 
k=\ 

as can be easily seen from the properties of the Fourier 
transforms. However, even if l(u) is not a polynomial, 
equation (fT^ defines an operator which is called pseudo- 
differential. We will now analyze the properties and 
the role of our pseudodiffcrential operator r]{dx) to see 
if (fTO|) can reasonably be considered as a generalized 
Schrodinger equation. 
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III. SEMIGROUPS AND GENERATORS 

Let X{t) be a one dimensional Levy process, namely 
a process with stationary and independent increments, 
and X{0) = almost surely. The chf of its increments 
on a time interval At then is [ip{u)]^*^'^ where f{u) is an 
infinitely divisible law, and r a time scale parameter (see 
for example [8j] and jUj for details about Levy processes). 
It is well known that t](u) — In ip{u) is the Ich of an 
infinitely divisible law if and only if it satisfies the Levy- 
Khintchin formula JJJ 



77(1*) = iju - 



^ Jr 



- 1 — iuxl[_i i]{x)'\ v{dx) 
(13) 

where 7, € R, is the indicator 0-1 function of the set 
A, and v{-) is the Levy measure of our infinitely divisible 
law, namely a measure on R such that i'({0}) = and 



{x^ A 1) v{dx) < +00 . 



The integrals involving the Levy measure v should then 
in general be calculated on R — {0} to take into account 
its behavior near y = 0. The triplet (7, (3, v) completely 
determines the Levy process and is also called its charac- 
teristic triplet. There are^ a few equivalent formulations 
of this important result [S] . In particular the truncation 
function -^[-1,1] (x) can be chosen in several different ways; 
this choice, however, will affect only the value of 7, while 
P and ly would be left unchanged. 

To every Levy process is associated a semigroup 
{Tt)t>() acting on the space T) of the measurable, 
bounded functions [TT'|: if / G P, we have (T*/) {x) = 
Elf {X{t) + x)], where E is the expectation. The in- 
finitesimal generator A of the semigroup (see fn^ p. 131) 
is now defined on the domain Va of the functions / G I? 
such that the limit (in norm on V) 



Af = lim 



Ttf-f 
t 



exists. It can be proved (see [ll| p. 139) that the gen- 
erators of a Levy process are pseudodifferential oper- 
ators that can be extended to the Schwartz space S 
of the rapidly decreasing functions. In particular we 
find that the symbol of A is nothing else than //(u), 
namely A ~ r]{dx), and that from the Levy-Khintchin 
formula we have 

[Vid.)f]{x)^j{d,f){x) + ^{dlf){x) 

+ I [f{^ + y)-f{x)^y{d.f){x)i[-i,i]{y)\ Hdy)- 



In other words our pseudodifferential operator rj{dx) of 
equation (jlOp is the generator of the underlying Levy pro- 
cess, and thanks to the Levy-Khintchin formula it also 
has an explicit expression in terms of integro-differential 
operators. 



The generator A = ri{dx) can finally be extended to 
i^(R) which is a Hilbert space, so that we can also dis- 
cuss its self-adjointness. In particular if X{t) is a Levy 
process, then its infinitesimal generator A = i]{dy) will 
be self-adjoint in L^(R) if and only if X{t) is centered 
and symmetric, namely if the symbol ri{u) is real with 



= — -u 



(cos ux — 1) ^{dx) 



(14) 



where ) is a symmetric Levy measure (see [ll| p. 154). 
A Levy measure is symmetric when i'{B) = v{—B) for 
every Borel measurable set, where —B = {x]—x e _B}. 
As a consequence the self-adjoint generators of the cen- 
tered and symmetric Levy processes enjoy the following 
simplified, integro-differential form 



{Af){x) = [ii{dx)f]{x) 



(15) 



{dlf){x)+ / [f{x + y)~ f{x)]v{dy). 

It is also possible to show that ~ri{dx) is not only self- 
adjoint, but also positive on L^(R) in the sense that for 
every / S L^(R) we have — (/, ?7(<9a;)/) > 0, and this is 
equivalent to say that the spectrum of —ri{dx) lies entirely 
in [0, -l-oo). 

We come back now to our problem of giving a general- 
ized, Levy-Schrodinger equation. Let X(t) be a centered, 
symmetric Levy process with pdf and transition pdf: we 
know that there is a centered, infinitely divisible law with 
chf ip(u) = e^i"^) such that the chf of its stationary in- 
crements AX{t) ^ X{t + At) - X{t) is 



At/r _ g?7(u)At/T 



for a suitable time scale parameter r. The transition pdf 
then is © , so that by means of the substitution t^itwe 
get the propagator G of equation ^ . A wave function ij} 
then evolves following (|10p . and since our process is cen- 
tered and symmetric the generator ri{dx) is self-adjoint 
and has the integro-differential expression (|15p . That 
means that our proposed equation takes the form 

idt^{x,t) = -"^^Pix^t) (16) 

T 

If 

= -^dliPix) - - j^i'{x + y)- i'ix)] v{dy). 

When X{t) is a Gaussian Wiener process we know that 

rj{u) = ~fP-u^l2 and A = r^ipx) = Then the 

process evolution equation is reduced to the Fokker- 
Planck equation (|4]) with D = /3^/2t, and we have ar- 
gued that in this case the usual non relativistic, free 
Schrodinger equation can de obtained by means of the 
substitution ([7]). In fact this amounts to take the - self- 
adjoint and positive in L^(R) - (pseudo)differential op- 
erator —hri{dx)/T — —hDd^ as the kinetic energy opera- 
tor. We propose here to extend this shortcut also to the 
case of non Gaussian Levy processes. Apparently this is 
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just a formal analogy, but there are at least two ways 
to make it more compelling: the Feynman path integral, 
and the Nelson stochastic mechanics. Here we only take 
for granted this association between generators and ki- 
netic energy operators in order to establish the form of 
the free Levy-Schrodinger equation and discuss its first 
consequences; we instead postpone to a subsequent pa- 
per a more rigorous derivation based on the application of 
Nelson stochastic mechanics to Levy processes. We will 
remember however that the method of Feynman path in- 
tegrals has already been used in the particular case of the 
fractional Schrodinger equations [20| to obtain the same 
association, albeit in a more restricted case: that of the 
stable laws, a particular class of infinitely divisible laws. 
Note that we will adopt here only the formal substitution 
t <-> it, but we will not impose D = h/2m because our 
generalizations are not necessarily supposed to be some 
kind of quantum mechanics, but will rather describe a dy- 
namical theory of Levy processes in the spirit of Nelson 
stochastic mechanics. We can also introduce a constant 
a with the dimensions of an action, so that our proposed 
free Levy-Schrodinger equation (jl6p becomes 



iadt-ipix, t) ^ Ha^p{x, t) = ri{dx)i'{x, t) 

T 



(17) 



-o^TT ^xV'(a;, t)--[ [ip{x + y,t)- ^{x, t)] v{dy) 

where now the free hamiltonian operator Hq has the di- 
mensions of an energy. This integro-differential hamilto- 
nian Hq is self-adjoint and positive on L^(R) so that it is 
a good kinetic energy operator. Everything that can be 
deduced about the usual Schrodinger equation from the 
positivity and self-adjointness of Hq can also be of course 
extended to our Levy-Schrodinger equation ^T7]\ . In par- 
ticular the conservation of the probability in the sense 
that, if IV'P plays the role of the position pdf, then the 
norm HV'P will be constant. We could finally add a po- 
tential V{x) to P?)) and get a complete Levy-Schrodinger 
equation 

iadt'ipixjt) = Hil){x,t) — 

T^{d.,)^P{x,t) + V{x)ij{x,t) (18) 



where the hamiltonian is now H = Hq + V. 



IV. DISCUSSION AND EXAMPLES 

The free Levy-Schrodinger hamiltonian of equa- 
tion (jl7p contains two parts: the usual kinetic energy 
— related to the Gaussian part of the process; and 

the jump part which is given by means of an integral 
with a symmetric Levy measure v. Of course, depend- 
ing on the nature of the underlying process, the equa- 
tion (|17p can contain these components in different mix- 
tures. If the underlying process is purely Gaussian then 
the Levy measure ly vanishes and (jl7p is reduced to the 



usual Schrodinger equation. On the other hand when the 
initial process is totally non Gaussian, then /3 = and we 
get a pure jump Levy-Schrodinger equation. In general 
both terms are present and if for instance we introduce 
Lu — 1/t and choose a = h and = ar/m, then pT|) 
takes the form 



ihdt'4>{x, t) 



2^ dx'>P{x, t)-hw J [ipix + y,t)- ipix, t)] v{dy) . 

Here the jump term can be considered as a correction to 
the usual Schrodinger equation and its weight, the energy 
Huj, is at present a free parameter. In particular for P — 
we get pure jump Schrodinger equations of the form 

idt^ix, t)^-uj [ [^(x + y,t)- ^(x, t)] u{dy) . (19) 
Js. 

Of course these remarks emphasize the fact that the ex- 
plicit form of the Levy-Schrodinger equation will depend 
on the choice of the particular Levy measure u. 



A. Stationary free solutions 

Let us consider first the stationary solutions of \n\ : 
taking 



i^{x,t) 



_ ^-lEt/a 



(j){x) 



iadtipixjt) = E^p(x,t) 



we have that the spatial part 4'{x) will be solution of 

Hocb{x) = -^<t>"{x)-- [[<t>{x + y)-^{x)]v{dy) 



= E(l){x) 
For the plane wave solutions 
(p{x) = : 



(20) 



and because of the symmetry of the Levy measure v, 
equation ([20]) becomes 



E(t>{x) 



(cos uy — 1) h'idy) 



--7]{u)(t){x) 
T 



and hence it is satisfied when between E and u the fol- 
lowing relation holds 

a 



E 



-■q{u) 



Here w is a wave number, while we are used to look for a 
relation between energy E and momentum p. If then we 
posit p = au the energy-momentum relation for our free 
Levy-Schrodinger equation is 



E = 



T \a/ 



(21) 
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B. Some classes of infinitely divisible laws 

We will explore now a few examples of Levy- 
Schrodinger equations associated to Levy processes. For 
a short summary of the concepts used here see for ex- 
ample [lOl- We will consider the c/i/'s, Levy measures 
and infinitesimal generators of centered, symmetric, in- 
finitely divisible laws so that p4p and (fT5|) hold. The 
form of the Levy measure v is then instrumental to ex- 
plicitly show how the pseudo-differential generator rjidx) 
works. It would then be useful to list several classes of 
infinitely divisible laws in a growing order of generality: 



1. Stable laws: here we have [sl. [l0| 
{a\u\)^ . 



7]{u) = -- 



A 



< A < 2, 



(22) 



with the important particular cases 

-a^v? /2 Gauss law (A = 2); 



77(u) 



-a\u\ 



Cauchy law (A = 1) 



Stable laws are selfdecomposable and hence their 
Levy measures are absolutely continuous [1, so 
that v{dx) — W{x) dx with 



W{x) 



B 



A+l 



B = 




(23) 



A = 2; 
A=l; 
A 7^ 1,2. 



The infinitesimal generator, which in the Gauss 
case (A = 2) simply is 

y(al/)(a:), 

for < A < 2 becomes the pseudo-differential op- 
erator 



I 



+ 00 



B 



2tT j^oo 

+00 



{a\u\y 
A 



e'''''f{u)du 



fix + y)- fix) 



\y 



A+l 



dy 



which can also be expressed symbolically in terms 
of the fractional derivatives [l9| 



1 



/27r J-oo 

^ f + OO 



|M|^e*"^7(w) du 
i-dlfl^e^--fiu)du 



\/27r J- 

The Cauchy case is discussed at length in [2l| 



2. Selfdecomposable (nan stable) laws: they are an 
important sub-family of infinitely divisible laws. 
Two ex amp les are (for details and other exam- 
ples see [23^'30']) the Variance Gamma family (VG, 
A > 0) and a law connected to the relativistic quan- 
tum mechanics: 



riiu) = 



-Al n(l + a^u^ ), VG; 

1 — Vl + a^u^ ^ relativistic q.m. 



which have no Gaussian part (/? = in the Lcvy- 
Khintchin formula) and produce pure jump pro- 
cesses. Their Levy measures have densities 0, l30l | 



Wix) 



A|a;|-le-l^l/^ VG; 
(7r|2;|)~^iiri(|a::|/a), relativistic q.m. 



where K\iz) is a modified Bessel function. Remark 
that, while for the VG pdf can be explicitly given 



qix) 



^2^r(A)\/27r V a 



we have no elementary expressions for the Rela- 
tivistic q.m. pdf 



qix) 



1 

2^ 



+ 00 



'du 



which can then be calculated only numerically. 

3. Infinitely divisible (non selfdecomposable) laws: 
The classical example of an infinitely divisible, non 
selfdecomposable law is the Poisson law of inten- 
sity A, but the corresponding Levy process would 
not be symmetric. If however we take the chf x(^) 
of a centered, symmetric law, then the correspond- 
ing compound Poisson process will be centered and 
symmetric with Ich 

riiu) = X[x{u) - 1] . 

In the analysis of the corresponding Levy measure 
we must remember that now we can no longer take 
for granted that v is absolutely continuous. For ex- 
ample if the jump size can take only two values ±a 
(a > 0) with equal probabilities 1/2, then the chf 
x(u) = cosau has no pdf rjiu) = A(cosau— 1), and 
vidx) — XFidx) where the cumulative distribution 



Fix) 



eix -a) + eix + a) 



(24) 



is a symmetric, two-steps function, and 0(a;) is the 
0-1 Heaviside function. If on the other hand x is a 
law with a pdf gix), it is possible to show that also 
the Levy measure v is absolutely continuous with 
a density 

Wix) = Xgix) . 

This completely specify the associated Levy process 
on the basis of the Poisson intensity A, and of the 
law of the jump sizes. 
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C. Examples of Levy— Schrodinger equations 

We can now analyze a few examples of Levy- 
Schrodinger equations based on the laws listed above. 
The first two cases have an apparent physical meaning, 
while at present the other three (Variance-gamma, stable 
and compound Poisson laws) are short of an immediate 
interpretation. 

1. Non relativistic, quantum, free particle: this is the 
well known case of the Gaussian Wiener process 
with 7/(m) = — /3^u^/2 giving rise to the usual 
Schrodinger equation ^ for a suitable identifica- 
tion of the parameters. In this case the energy- 
momentum relation (|2ip is 



[3^ 



0^ 
2aT 



P 



and with a = h and = olt jm we get as usual 

r,2 



E = 



V 
2m 



p — hu . 



2. Relativistic, quantum, free particle: it is interesting 
to remark at this point that there is a Levy process 
which is connected to the relativistic Schrodinger 
equation, in the same way as the Wiener process is 
connected to the non relativistic Schrodinger equa- 
tion. Take the non stable, selfdecomposable law 
ri{u) = 1 — Vl + a^?!^, and use the following iden- 
tifications 



a 2 
— = mc 

T 



— , p = hu 

mc 



to find from (EI 



E 



-mc^ri{u) — \/ m^c'^ -1- p^c^ — m(? (25) 



which is the relativistic total energy less the rest 
energy mc?: namely the kinetic energy. The corre- 
sponding Levy-Schrodinger equation is now 



ihdti^{x,t) — \/ nv^d^ — c^f?d'^ — m(? ip{x,t) 

and is discussed at length in ^2l|. Since the con- 
stant —mc^ can be reabsorbed by means of a phase 
factor e*™^ the wave equation finally is 



ihdtip{x,t) — ^m'^c^ — (?-W-d'^i\;{x,t) 



(26) 



which is the simplest form of a relativistic, free 
Schrodinger equation ^33*1. It is interesting to note 
that the Levy process behind the relativistic equa- 
tion (j26[) is a pure jump process with an absolutely 
continuous Levy measure with pdf 



Tr\x\ \ / ""fI V h 



so that equation ((26)) can also be written as 
ihdt^{x,t)= (27) 

From the form of the relativistic energy (PS)) also 
the usual relativistic corrections to the classical 
energy-momentum relation for small values oi p/c 
follow: 



E 




p- 



P 



o{p' 



2m Sm^c"^ 

Finally if we consider E — H{p) as the Hamilto- 
nian function of a relativistic free particle from the 
Hamilton equations we get 

q = OpH = 

m -^p'^jm^c^ 

and here too we can see the relativistic correction 
to the classical kinematic relation p — mq. 



3. Variance- Gamma laws: for A 
1 



i we have 



r]{u) 



-ln(l + aV), g(x)=.-i^o(^ 
2 an \ a 



W{x) 



'\x\/a 



2\x\ 



The Variance-Gamma processes are pure jump pro- 
cesses with no Gaussian part in the Levy-Khintchin 
formula (fT3|) {/3 — 0), so that the Levy-Schrodinger 
equation becomes 

iadt'ip{x, t) — 

_Xa r i;{x + y,t)-^{x,t) ^_|^|/,^^ ^28) 

r Jr \y\ 

By choosing 



ma 
p = u 



' r (29) 

T r 

we have the following energy-momentum relation 
(for pr/ma 0) 



E 



2t2 
P' 



In 1 



T^p^ 



t2 p4 



+ o(/) 



while with the identification E — H{p) we can also 
recover the kinematic relations between p and q: 



q 



p/m 



P 



1 



^2p2 



o(/) 



It is apparent then that again these equations are 
corrections to the classical relations. 
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4. Stable laws: for the non Gaussian stable laws we 
have the Ich (|22p and the Levy measure pdf ([^5]) 
where < A < 2. With a couple of exception 
(Cauchy and Levy laws), however, there are no ele- 
mentary formulas for their prf/'s q{x). Their Levy- 
Schrodinger equation is 



iadtip{x, t) 



T 2Acos^r(-A) 



\y 



A+l 



(30) 
dy 



the energy- 



and with the identifications 
momentum relations become 

,„2\ 1-^/2 / 9 X A/2 



T A 



2A/2 



2m J 



This however looks not as a correction of the clas- 
sical formula as in the other cases considered, but 
rather as a completely different formula. The same 
can be said of the kinematic relations between q 
and p which arc now 



P 

9 = — 



2 2 \ A/2-1 



These relations are still another reason to consider 
not advisable to restrict an inquiry on Levy pro- 
cesses and Schrodinger equation only to the family 
of stable processes. 

5. Compound Poisson process: for the symmetric, 
compound Poisson process with the Levy mea- 
sure given by (|24p the pure jump Levy-Schrodinger 
equation (fTO|) greatly simplifies as 

idt^p{x,t) = 

— ^ ['/'(a; + a, i) - 2ip{x, t) + ip{x - a, t)] . 

We can now show that the usual Schrodinger equa- 
tion ([T]) can always be recovered as a limit case of 
this Poisson-Schrodinger equation. If ip is twice 
differentiable in a;, we know that for a ^ 0+ 

ijjix ± a, i) = 

^{x, t) ± ad^)i^{x, t) + —dlij{x, t) + o(a2) 
and hence 
idttp(x,t) 



dl'4>{x,t)+\o{a^). 



Now, if, as a — )■ 0^, also A +oo in such a way 
that Aa^ — > fe^, then we have 

ujb"^ 

idMx.t) = — —dl^{x,t) 

where w6^/2 has the dimensions of a diffusion co- 
efficient: namely, in the limit, we get a Wiener- 
Schrodinger equation of the type (P). This proce- 
dure can also be used to introduce small corrections 
in the coefficient It? /2m of a Schrodinger equation. 



D. Perfectly rigid walls 

An example of solution of the complete Levy- 
Schrodinger equation (fTS]) can easily be obtained in the 
case of a system confined between two perfectly rigid 
walls symmetrically located at a; = ±i/2. The discussion 
is similar to that of Section flV Al but for the boundary 
conditions which now require that the solutions vanish at 
X ^ ±L/2. As a consequence the solutions are the usual 
trigonometric functions with discrete eigenvalues 



En 



'—V{Ur. 
T 



mr 

— , n^l,2, 



The form of the eigenvalue sequence will depend on the 
Ich rj. If the underlying process is a Wiener process we 
have ri{u) — — /3^u^/2, and hence 



En — 



which, with the identifications a — h and = ar/m, 
coincides with the usual quantum mechanical result. On 
the other hand for an underlying Variance-Gamma noise 
we find 



In 1 



Er, 



Finally for a symmetric, compound Poisson noise with 
intensity A and equiprobable jump sizes ±a we get 



a7rn\ 



En = 1 - COS ■ 

r V L / 

Apparently this no longer is a monotone increasing se- 
quence: En goes up and down between and 2aA/T. If 
a/ L is rational the sequence is periodic; on the other hand 
when a/ L is irrational there are no coincident eigenvalues 
in the sequence, so that the i?„ will fill the bandwidth 
between and 2aA/r. 



V. CONCLUSIONS 

We have discussed the possibility of generalizing the 
relation between Brownian motion and Schrodinger equa- 
tion by formally associating the kinetic energy of a more 
general system to the generator of a symmetric Levy pro- 
cess, namely to a pseudodifferential operator whose sym- 
bol is the Ich of an infinitely divisible law. This amounts 
to suppose, then, that this new Levy-Schrodinger equa- 
tion is based on an underlying Levy process that can have 
both Gaussian and jumping components. 

In recent years other extensions of the Schrodinger 
equation have been put forward in the same spirit of our 
Levy-Schrodinger equation. In particular we refer to sev- 
eral papers about fractional Schrodinger equations (20| 
that explored the use of fractional calculus in a kind 
of generalized quantum mechanics. As it is clear from 
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the previous sections, however, this is the particular case 
when our underlying process is stable. The extension to 
the infinitely divisible, non stable processes, on the other 
hand, is meaningful because there are significant cases 
that are now in the domain of our Levy-Schrodinger 
picture. In particular, as shown also in [2ll |. the sim- 
plest form of a relativistic, free Schrodinger equation can 
be deduced from a particular type of selfdecomposable, 
non stable process. Moreover in many instances the new 
energy-momentum relations can be seen as corrections 
to the classical relations for small values of certain pa- 
rameters. It must also be remembered that our model 
is not tied to the use of processes with infinite variance: 
the variances can be chosen to be finite even in a purely 
non Gaussian model - as in the case of the relativistic, 
free Schrodinger equation ~ and can then be used as a 
measure of the dispersion. 

It is important now to explicitly give in full detail a 
derivation of the Levy-Schrodinger equation from either 
the Feynman integrals or a generalized stochastic me- 
chanics. This seems to be possible because the techniques 
of the stochastic calculus applied to Levy processes are 
today in full development ^ [l^, [ll|, Il2l |. and at our 
knowledge there is no apparent, fundamental impediment 
along this road. At present our approach lacks this rig- 



orous discussion of how the Levy-Schrodinger equation 
comes out from the evolution equations of the Levy pro- 
cesses. We have confined ourselves to give only a few 
heuristic arguments based on both the identification of 
the process generators as the kinetic energy operators, 
and the analytic continuation of the time variable t to 
its imaginary counterpart it. We hinted, however, to two 
possible ways of giving a more rigorous derivation: we 
can first of all follow the Feynman integral road. In this 
case we should bear in mind that the relations among 
kinetic energy and momentum are no longer the usual 
relations: this is important to correctly write the La- 
grangian in the Feynman integral. Alternatively we can 
try to generalize Nelson stochastic mechanics by adding 
a suitable dynamics to our Levy processes, and this will 
be the subject of a future paper. 
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